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^ ! Abstract 

This paper presents a method for approximating a class of complex transfer func- 
^ ■ tion matrices corresponding to physically realizable complex linear quantum sys- 

C/^ ■ tems. The class of linear quantum systems under consideration includes intercon- 

c/3 ! nections of passive optical components such as cavities, beam-splitters, phase- 

shifters and interferometers. This approximation method builds on a previous 
result for cascade realization and gives good approximations at low frequencies. 

J:!^ ■ Keywords: Quantum Linear Systems, Model Reduction, Cascade Realization. 
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. . 1. Introduction 

Q . In recent years, there has been considerable interest in the modeling and feed- 

back control of linear quantum systems; e.g., see jl MlTll . Such linear quantum 



systems commonly arise in the area of quantum optics; e.g., see [|l8l Il9l1 . The 



L^ ■ feedback control of quantum optical systems has applications in areas such as 

'^ ■ quantum communicatio ns, quan tum teleportation, and gravity wave detection. In 



9^. , particular, the papers ISi llSMlTII have been concerned with a class of linear quan 

tum systems in which the system can be defined in terms of a set of linear complex 
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quantum stochastic differential equations (QSDEs) defined purely in terms of an- 
nihilation operators. Such linear complex quantum systems correspond to optical 
systems made up of passive optical components such as optical cavities, beam- 
splitters, and phase shifters. This paper is concerned with the approximation of 
systems in this class of linear complex quantum systems. The method proposed in 
this paper builds on the result of ISiilTl] which gives a method for physically real- 
izing a given complex transfer function matrix corresponding to a linear quantum 



system in the class considered in rtl5L ll6ll. 



In the approximation of linear quantum systems, it is important that the ap 
proximate system which is obtained is physically realizable. The issue of physi 



cal realizability for linear quantum systems was considered in the papers 141151.110. 



15L116D. This notion relates to whether a given QSDE model represents a physical 
quantum system which obeys the laws of quantum mechanics. In applying apply- 
ing approximation methods to obtain approximate models of quantum systems, it 
is important that the approximate model obtained is a physically realizable quan- 
tum system so that it exhibits the features inherent to quantum mechanics such as 
the Heisenberg uncertainty principle. 

The approximation method proposed in this paper follows directly from the 
physical relation algorithm proposed in [9l, llTO. The physical realizability of the 
approximate system follows directly from the fact that the algorithm proposed in 
[Isi [iTD leads to a physical realization in terms of a cascade connection of optical 



cavities. For this approximation method, we present some bounds and approxi- 
mate bounds on the approximation error as a function of frequency. 

One application of the approximation method proposed in this paper is in mod- 
elling of linear quantum systems where it is desired to construct a simpler, but still 
physically realizable model of a complex quantum linear system in such a way that 
a frequency dependent bound on the approximation error is obtained at low fre- 
quencies. Another application of the approximation method occurs in the case of 
coherent quantum feedback control systems when both the plant and controller are 
linear quantum systems; see [lli|2i|4i[la|20|]. In this case, it is desired to construct 
a simpler coherent quantum controller which is still physically realizable. 

2. A Class of Linear Complex Quantum Systems 

We consider a class of linear quantum systems described in terms of the anni- 
hilation operator by the quantum stochastic differential equations (QSDEs): 

da{t) = Fa{t)dt + Gdu{t); 

dy{t) = Ha{t)dt + Jdu{t) (1) 



where F e C"^", G 
Hliillil]. Herea(t) 



Oi (t) ■ ■ ■ a„(t)] is a vector of (linear combinations of) 
annihilation operators. The vector ^(t) represents the input signals and is assumed 
to admit the decomposition: du{t) = f3u{t)dt + du(t) where u{t) is the noise part 
of M(t) and Pu{t) is an adapted process (see UM, ll23ll and nlm ). The noise u{t) is 
a vector of quantum noises. The noise processes can be represented as operators 
on an appropriate Fock space (for more details see ll25ll . t2M)- The process (3u{t) 
represents variables of other systems which may be passed to the system ([T]) via 
an interaction. More details on this class of systems can be found in [.15.1 . L4J. 

Definition 1. (See M lUu-j A linear quantum system of the form ([7]) is said to be 
physically realizable if there exists a commutation matrix = 0^ > 0, a coupling 
matrix A, a Hamiltonian matrix M = M^, and a scattering matrix S such that 



-0 



iM + ^AU j ; G 



-eA^S-H = A-J = S 



(2) 



and S'^S = I. 



Here, the notation ^ represents complex conjugate transpose. In this definition, 
if the system ([U) is physically realizable, then the matrices S, M and A define 
a complex open harmonic oscillator with scattering matrix 5", coupling operator 
L = Aa and a Hamiltonian operator "H = a^Ma; e.g., see pill. 112311 . [|22D, BD and 



i,[l5 



1511 to allow for 



12611 . This definition is an extension of the definition given in 
a general scattering matrix S; e.g., see [J6l,ll2i]. 

The following theorem is an straightforward extension of Theorem 5.1 of lll5n 
to allow for a general scattering matrix S. 

Theorem 1. (See ^13].) A complex linear quantum system of the form ([7]) is phys- 
ically realizable if and only if there exists a matrix = 0^ > such that 



FQ + QF^ + GG^ = 0; G = -QH^J; JU = L 
In this case, the corresponding Hamiltonian matrix M is given by 



(3) 



and the corresponding coupling matrix A is given by 

A = H. 



(4) 



(5) 



Definition 2. The linear complex quantum system ([7]) is said to be lossless bounded 
real if the following conditions hold: 

i) F is a Hurwitz matrix. 

ii) $(s) = H{sl - F)-^G + J satisfies <l>(iw)t<l>(iw) = I for all w G M. 

The following definition extends the standard linear systems notion of minimal 



realization to linear complex quantum systems of the form ([T]); see also [llSh . 

Definition 3. A linear complex quantum system of the form ([7]) is said to be min- 
imal if the following conditions hold: 

i) Controllability. a^F = \a^ for some A G C and a^G = implies a = 0; 
ii) Observability. Fa = \afor some A G C and Ha = implies a = 0. 



The following Theorem is an straightforward extension of Theorem 6.6 of ||15|] 
to allow for a general scattering matrix 5*. 

Theorem 2. A minimal linear complex quantum system of the form ([7]) is physi- 
cally realizable if and only if the system is lossless bounded real. 

Definition 4. The complex linear quantum system ([7]) is said to be a quantum 
system realization of a complex transfer function matrix K{s) if 

K{s)=H{sI-F)-^G + J. (6) 

3. The Cascade Cavity Realization Algorithm 

In this section, we recall the cascade cavity realization result of ^ ilTD and 
generalize it slightly to allow for quantum systems with a more general scatter- 
ing matrix. Indeed, given a linear quantum system of the form ([T]) with transfer 
function matrix ©, we can write K(s) = JK{s) where 

K{s)=H{sI-F)-^G + I (7) 

and H = J^^H. Corresponding to © is the linear quantum system 

da{t) = Fa{t)dt + Gdu{t)] 

dyit) = Ha{t)dt + du{t). (8) 



In order to obtain a physical realization of ^, the result of |19ill7|] can be applied 
to transfer function matrix (|7]). Then a collection of beamsplitters can be used to 
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Figure 1: Cascade of n, generalized m mirror cavities and a unitary operator. 



implement the unitary matrix J; e.g., see [27]. This leads to a physical realization 
of the transfer function matrix (|7]) as shown in Figure [TJ 

An optical ring cavity consists of a number of partially reflecting mirrors ar- 
ranged to produce a traveling light wave when coupled to a coherent light source; 
e.g., see 11191. 12111. If we augment such a cavity by introducing phase-shifters on 
the input and output channels, such a cavity with m mirrors, can be described by 
a linear quantum system of the form ([8]) as follows; see [191. 1 1711: 



da 



padt — h^du;dy = hadt + du 



where 



p+p* = -1 = -Yl 



Ki 



-h^h. 



(9) 



(10) 



2=1 



Here p = —7/2 + iA, 



h 



hi 
h2 



Kie 



K2e 



rirn,^ 



J6»2 



, du 



dui 
du2 

GjUim. 



dy 



dyi 
dy2 

dym 



Furthermore, any first order complex linear quantum system of the form ^, with 
non-zero h E C™ and satisfying (flOl) . can be physically realized as a generalized 
m mirror cavity. In this case, the mirror coupling coefficients and phase shifts are 
determined using a polar coordinates description of the elements of h. Also, the 
detuning parameter A is determined from the imaginary part of the system pole p. 
The cascade cavity realization introduced in [|9|, [l7|] involves a cascade in- 
terconnection of n, generalized m mirror cavities as shown in Figure [2l In this 
cascade system, the ith cavity is described by the following QSDEs of the form 
I®,®: 



dtti 



Pittidt — Hjdu; dy = Hiaidt + du 



(11) 
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Figure 2: Cascade of n, generalized m mirror cavities. 



where 



Pi+P^ 



-H]m. 



(12) 



The cascade system is then described by a complex linear quantum system of the 
form ([8]) where 



Vi 







H 



[H, H2 ... Hn],J = I. 



G = - 



Hl 
Hi 



Hi 



(13) 



Reference [J9i LLTI] presents an algorithm for realizing a physically realizable 
quantum system ([8]) with transfer function (|7]) via a cascade of generalized cavi- 
ties. We restrict attention to quantum systems in which the transfer function (|7]) 
corresponds to a minimal system ^ such that the eigenvalues of the matrix F are 
all distinct. In this case, it follows via a (complex version of a) standard result 
from linear systems theory, that the system ([8]) can be transformed into Modal 
Canonical Form; e.g., see uM . The complex linear quantum system in modal 
canonical form is assumed to be as follows: 



where 



da{t) 


= Fa{t)dt + 


Gduit) 


dy{t) 


= Ha{t)dt + du{t) 


Pi ... ■ 




[Gi 1 




P2 ■ 


;G = 


G2 


;H = 


': ■•. 








... pn. 




[gJ 





(14) 



Hi H2 : Hn 



(15) 



Also, it is assumed that in this realization, the eigenvalues are ordered so that 

bi| < b2| < • • • < |Pn|- Then, K{s) satisfies the equation 



K{s) = H{sI-F)-^G + I. 
The algorithm proposed in [|9|,[l7|] is as follows: 



(16) 



Step 1: Begin with a minimal modal canonical form realization (fT4l) . (fTSJ) of the 
lossless bounded real transfer function matrix K(s). 



Step 2: Let 



H„ 



H„ 



H a 



Pn+Pn 



H„ 



lOLr. 



t{n,n) 



'Oir. 



(17) 



Step 3: Calculate the quantities Hn, Hn-i, ■ ■ ■ , Hi, an, Oin-i, . . . , ai, t{i,j), for 
j = n,n — 1, . . . ,1 and j > i. These are calculated using the following 
recursive formulas starting with the values determined in Step 2 for i = n: 



Hi 



Oli 



^+E 



Hi 



j=i+ 



iP^-P^t 



E ^(^■' ^)^i 



--i+i 



Hi\ 



.Iilk 

Pi + pI 



-, H^ 



H 



i^i 



t{k,i) 



t{i,i) 



1 ^ 

- — J2 i{k,])H]Hi for fc = 2 + 1, . . . 



Pi -Pk 



j=i+i 



'OLi 



(18) 
(19) 



(20) 
(21) 



Step 4: Set t{k, i) = for k < i and define an n x n transformation matrix T 
whose (i, j)th element is t{i,j). 

The following theorem is presented in liOlllTn. 



Theorem 3. Consider anmxm lossless bounded real complex transfer function 
matrix K{s) with a minimal modal canonical form quantum realization 47?1) . ( |75l) 
such that the eigenvalues of the matrix F are all distinct and that all of the matrix 
inverses exist in equation rfTSl) when the above algorithm is applied to the sys- 
tem f l7?l) . d75l) . Then the vectors Hi, H2, . . . , Hn defined in the above algorithm 
together with the eigenvalues pi, p2, . . . ,Pn define an equivalent cascade quan- 
tum realization ([S]), ((Til) for the transfer function matrix K{s). Furthermore, this 
system is such that the condition 472]) is satisfied for all i. Moreover, the matrices 
{F,G,H,I} defining this cascade quantum realization are related to the matrices 
{F, G, H, 1} defining the modal quantum realization d7?l) . f lTTl) according to the 
formulas: 

F = TFT-\ G = TG, H = HT-^ (22) 

where the matrix T is defined in the above algorithm. 

The physical realization of Q corresponds to writing 

K{s) = Jkn{s)kn-i{s) . ..ki{s) (23) 

where each transfer function matrix Ki{s) is a first order transfer function matrix 
corresponding to an optical cavity described by a QSDE of the form (fTTI) . 

4. The Main Result 

Our proposed method for obtaining an approximate model for a complex lin- 
ear quantum system ([T]) with transfer function matrix Q involves truncating the 
cascade realization (l23l) to obtain the approximate transfer function matrix 

Ka{s) = Jakr{s)kr-l{s) . . . ki{s) (24) 

where Ja = JKn{0)Kn-i{0) ■ ■ ■ -^r+i(0) and r < n is the order of the approx- 
imate model. It follows from this construction that Ka{s) is lossless bounded 
real and hence physically realizable. Indeed, since the transfer function ma- 
trix JKn{s)Kn^i{s) . . . Kr+i{s) is losslcss boundcd real, it follows that the ma- 
trix Ja will be unitary. Therefore, Kais) will be lossless bounded real since 
Kr{s)Kr-i{s) . . . Ki(s) is lossless bounded real. 

In order to construct a state space realization of the reduced dimension transfer 
function matrix Ka{s), note that it follows from the development in Section [3] that 



the transfer function matrix Ka{s) = Kr{s)Kr^i{s) . . . Ki{s) has a state space 
realization of the form ([8]) defined by the matrices 



Pi 







Hr, 



[Hi H2 ... Hr ] . 



; Ga 



Hi 
Hi 



Hi 



(25) 



Also, the transfer function matrix Kb{s) = Kn{s)Kn-i{s) . . . Kr+i{s) has a state 
space realization of the form dS]) defined by the matrices 



Pr+l 







—H,,.j^2Hr+l Pr+2 











m 



—HlHr+l . . . —H^Hn-l Pn 

[ Hj.^1 Hj.^2 ■ ■ ■ Hn \ . 



, Gb 



hUi 

hU 



Hi 



(26) 



Hence, the matrix Ja is given by Ja = -^ (-^ — H^F^ ^Gb) and the reduced dimen- 
sion transfer function matrix Ka{s) has a state space realization of the form ([T]) 

defined by the matrices (Fa,Ga, JaHa, Ja ) ■ 

The ordering of the eigenvalues in the cascade realization (fT5l) means that this 
model is expected to be a good approximation of original model at low frequencies 
CO « \pr+i\- The corresponding error system is defined by 



KJs) 



K{s)-Ka{s) 



J 



Kn{s)Kn-l{s) . . . kr+l{s) 

-ir„(0)ir„_i(0)...K,+i(0) 



k,.{s)K,.i{s)...ki{s){21) 



We now present a result which bounds the induced matrix norm of the approx- 
imation error \\K(,{jijj) \\ as a function of frequency. This bound will be defined in 



terms of the following quantities: 



iPi + Pi 



'^i^ibil-bi-j^l 



n n 



B2{uj) = u Y^ Y^ C,,^i,{ujy, 



ii=r+l J2=n+1 



Bkiu) 



n n 



Jl=r+1 J2=n + 1 «fc=ifc-i + l 



Bn^Auj) 



uC, 



r+l,r+2 



,...,n(^) 



(28) 



where 



Cji,J2,-,ifc('^) 



^ ■' ' '^m.^— 1 ^--^rn.2— TJ^l +1 ■'^7723 ^m2 + l '^^rrii ^i-iixo '^'^m'x 









nLi-(p.;+p*) 



nf^i bijMPi, -i'^l 



(29) 



Theorem 4. Consider a physically realizable linear complex quantum system of 
the form ([7]) and corresponding transfer function matrix ©. Suppose this system 
has a cascade cavity realization f [2?]) and a corresponding approximate transfer 
function matrix Ka{s) defined in ( [24\) . Then Ka{s) is physically realizable and 
the corresponding approximation error transfer function matrix Ke{s) defined in 
i{27\) satisfies the bound 



\KAjuj)\\<J2b. 



k{UJ 



(30) 



k=l 



for all oj > where the quantities Bk{(jj) are defined in i{28\) , (|29l). 
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Proof. The fact that Ka{s) is physically realizable follows from its definition 
as discussed above. Now for any w > 0, it follows from dTT] ) that 



\Ke{3^)\\ < \\J 



kn{jUj)Kn_i{ju) . . .Kr+l{iu) 

x\\kriju)kr^i{ju) . . .i^i(ja;)|| 

= II k4jUj)kn-lijiu) . . . kr+l{jUj) - kn{<^)kn-M . . . i^.+l(0) || 

(31) 

using the fact that the transfer function matrix kr{s)kj.-i{s) . . . ki{s) is lossless 
bounded real and the matrix J is unitary. 

We now consider the transfer function matrices 

k,{s) = kn{s)kn-l{s)...kr+l{s)-kn{0)kn-l{0)...kr+l{0) 

and ke{s) = kn{s)kn-i{s) . . . kr+i{s). Also, note that it follows from (fTTl) that 

"" — IT. fJ^ 

each transfer function matrix Ki{s) is of the form Ki(s) = I -\ ^. Hence, 



Pi-s 



we can write ke{s) = YVi=n ( ^ "I — 37 ) • From this it follows that we can write 



k,{s) = J2iUs)-Uo)) 



(32) 



fc=i 



where the transfer function matrices Ti(s), T2(s), . . . , Tn-r{s) are defined as fol- 
lows: 



Ti{s) 



n(s) 



Pi- s' 



E 

i=r+l 

il=r+l i2=*i+l 



{Ph - s) {pi 



il=r+l 12=11+1 ifc=ifc_i+l l=k ' 



Tn-r{s) = n^^^ 
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Now it follows from (|32l ) and the triangle inequality that 

n—r 

\\k,{juj)\\<Y,\\Tk{juj)-T,m 



k=l 



for all cj > 0. We now consider each of the terms Tk{juj) = Tk{juj) 
Indeed, for any w > 0, we obtain 



(33) 

Tfc(O). 



fiUu^) = E 



juH.HJ 



Pi (pi - M 

n 



i=r+l ■ 

n n 



n=r+l i2=n+l 



Tk{j(^) 



n n 



E E ••• E X{H^,Hl3^S..,....A3^): 

Jl=r+1 J2=«i+1 ik=ik-\+'^ l=k 



r+1 



Tn^rU^) = WH,H}iUjSr+l,,„.n{.3^) 



where 



Sn,i2,...,ikiJ^) 



(-ju;) +(-Jc^) 22m^iPim + 

/ ■ \k — 4 v^ k v-v k s~^ k n ■ n T) - 



( _q, ,'\^~P~^ ^^ spk s-^k T-rP 



/_\fc — p— Iv^fc sr^k yfc T-rfc — 1 

' -^ -^ -^m-i^— 1 ^m'2^'>Ti^-\-l'" ^'m'i._i — TTi^_2 + l -^^g — 1 "^^ 



nf=i P'; ('''i -•"") 



(34) 



(35) 
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Now it follows from ([HI) that ||-ffi-ff||| = HJHi = -pi - p* for i = r + l,r + 
2, . . . ,n. Hence, applying the triangle inequality and the Cauchy-Schwartz in- 
equality to (|34|) . it follows that for any cu > 

\\fi{juj)\\ < B.ioo); 
mUcon < B^iu); 

\\fk{jco)\\ < Bk{u)- 

||T„_,(ja;)|| < Br,-r{oo). (36) 

These bounds combined with (l33l) and (|3T1) lead to the inequality (|30l) . □ 



Remark 1. The quantity X]r=i Bk{uj) in bound 0O|) is probably too complicated 
to calculate for all cases except when n—r is equal to one or two. However, we can 
obtain some good approximations to this quantity which apply at low frequencies. 
Indeed, for u « \pr+i\, we obtain 

fc=l i=r+l l-P^I-l-P* ^^1 

for all w > 0. Furthermore since sup^>Q \pi — iu)\ = — |(pi + P*), we obtain the 
following useful upper bound on Bi{u): Bi{uj) < 2uj Y17=r+i T~\f'^^ '^^^ uj >Q. 

5. Illustrative Example 

To obtain our initial physically realizable quantum system, we start with a 
random Hamiltonian matrix Af > 0, a random coupling matrix A, and a random 
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commutation matrix 6 > defined as follows: 



M 



3.3314 

2.5448 - 0.8204J 
2.4007 -1. 1592j 
3.6470 -1.3066J 
1.9949 - 1.8970J 



2.5448 + 0.8204J 
3.3994 

2.7136 - 0.4185J 
3.7009 + 0.1246J 
2.6090 - 1.4039J 

3.6470 + 1.3066J 
3.7009 -0.1246J 
4.2612 - 0.3611J 
5.9568 
4.1173-1. 3450 j 



2.4007 +1.1592J 
2.7136 + 0.4185J 
4.1258 

4.2612 + 0.3611J 
3.5647 - 1.0224J 

1.9949 + 1.8970J 
2.6090 + 1.4039J 
3.5647 + 1.0224J 
4.1173 + 1.3450J 
3.9970 



e 



1.8356 

2.3408 + 0.3287J 
1.8732 - 0.0757J 
1.8750 - 0.0488J 
1.5306 - 0.3665i 



2.3408 
3.9779 
3.3007 
2.8626 
2.3995 



0.4982J 
0.2561i 
l.llGSi 

1.8750 + 0.0488J 
2.8626 + 0.2561J 
2.5536 - 0.7201J 
3.0974 
1.9480 - 0.9779J 



0.3287J 1.8732 + 0.0757^ 
3.3007 + 0.4982i 
3.8582 

2.5536 + 0.7201J 
3.1242 - 0.8781i 

1.5306 + 0.3665J 
2.3995 + 1.1 168j 
3.1242 + 0.8781J 
1.9480 + 0.9779J 
3.0319 



A = 



-1.0106 + O.OOOOj 0.5077 + 1.0950j 0.5913 + 0.4282J 
0.6145 - 0.3179i 1.6924 - 1.8740j -0.6436 + 0.8956J 

0.3803 + 0.7310J -0.0195 + 0.0403J " 
-1.0091 + 0.5779J -0.0482 + 0.6771J 

Also, we chose the scattering matrix S = I. This leads to a corresponding system 
of the form ([T]) where the matrices are defined as in (O. The eigenvalues of the 
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resulting matrix F are s = -0.0038 - O.OlSlj, s = -0.2103 - 0.1040J, s = 
-0.1674-1. 1066j, s = -3.0388-0.9275J, s = -10.8541-225.9473J. Clearly, 
the last eigenvalue has a much larger absolute value than all of the others and so we 
will apply our algorithm to approximate this fifth order system by a fourth order 
system. Bode plots comparing the original system frequency response with the 
reduced dimension system frequency response are shown in Figures [THU These 
Bode plots indicate that the proposed method gives a good approximation at low 
frequencies. Also, it follows from the construction that the reduced dimension 
system is lossless bounded real and so physically realizable. 



Bode Diagram 





10 10 

Frequency (rad/sec) 



Figure 3: Bode plot of original and approximate system from input 1 to output 1. 

In Figure |7l we show the singular value plot of the error transfer function 
matrix Ke{s) = K{s) — Ka{s) along with the error bound defined by Bi(uj). In 
this example, we see that the error bound is in fact exact since we only reduced 
the dimension of the original system by one. 

6. Conclusions 

In this paper, we have presented a method of approximating a class of linear 
complex quantum systems in such a way that the property of physical realizability 
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Bode Diagram 
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Figure 4: Bode plot of original and approximate system from input 2 to output 1 . 



Bode Diagram 
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Figure 5: Bode plot of original and approximate system from input 1 to output 2. 
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Bode Diagram 
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Figure 6: Bode plot of original and approximate system from input 2 to output 2. 



Singular Values 




* * * * 
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Figure 7: Singular value plot of the error transfer function matrix Kq{s) and the error bound 

Bi{uj). 
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(which is equivalent to the strict bounded real property in this case) is preserved. 
The paper presents a bound on the approximation error which shows that the ap- 
proximation is accurate at low frequencies. 
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